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Multiquark correlations inside hadrons can have a significant and in some cases even striking 
_C impact on the hadron spectrum. We show how such correlations in general, and mesons with a 

dominant tetraquark content in particular, emerge holographic ally in the AdS/QCD framework. 
On this basis, we arrive at a holographic realization of an exceptionally strong four-quark binding 
<^ and a correspondingly large tetraquark component in the lightest scalar mesons. Higher-lying 

tetraquark excitations, on the other hand, become too broad to form supernumeral scalar states. 
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1. Introduction 

A dominant tetraquark component is often invoked to explain non-standard features of mesons 
ranging from the lightest scalar nonet [1] to some of the recently discovered heavy-quark reso- 
nances [2]. Without guidance from the naive quark model, the theoretical analysis of such ex- 
otic states has to deal with the challenges of the strongly-coupled gauge dynamics more directly. 
Adding to the arsenal of suitable methods for this purpose (which now includes devoted lattice sim- 
ulations [3]), the AdS/QCD approach [4] was recently extended to provide holographic insights into 
tetraquark properties [5] . 

The AdS/QCD program describes strong-interaction physics analytically by constructing ap- 
proximate five-dimensional gravity duals for QCD on the basis of the gauge/string correspondence 
[6]. To access multiquark correlations inside hadrons within this framework, one exploits the fact 
that such correlations manifest themselves in the quark composition of the QCD interpolator which 
most strongly couples to the hadron under consideration. Hence these correlations leave a gauge- 
invariant imprint on the interpolator's anomalous dimension which the AdS/CFT dictionary [6] 
translates into a multiquark-content-dependent mass correction for the hadron's dual bulk mode. 
This holographic mechanism was previously employed to analyze the impact of diquark correla- 
tions on masses and sizes of the light-quark baryons [7]. In the following, we will apply it to 
tetraquarks. 

2. Tetraquark holography 

The AdS/CFT dictionary associates the bulk mode (p(x,z) (where x are the coordinates along 
the boundary and z denotes the orthogonal direction) dual to a given hadron \h) with the gauge- 
theory interpolator J(x) which most strongly couples to \h). More specifically, for a scalar meson 
the scaling dimension A of / determines the boundary condition [6] 



for the dual bulk-mode solutions. This condition is imposed on the (at small z leading) solutions of 
the bulk field equation by adjusting its mass term to 



(where m=, is the mode's mass and R is the AdSs curvature scalar). Interpolators with a larger quark 
content and correspondingly larger A are therefore related to heavier bulk excitations and generally 
to mesons with larger mass. (Note that the relation (2.2) is independent of the specific gravity dual 
under consideration). 

In contrast to the unique quark-antiquark interpolator J A = q a t A q a for ordinary scalar mesons 
with Aq q = 3, there are several independent (gauge-invariant and local) scalar four-quark oper- 
ators with a common classical scaling dimension A^z q 2 = 6. (An explicit example is J^ 2 = 
e abc e ade q h CT A q c q d CT A q e which contains a "good" (i.e. maximally attractive) diquark and a good 
antidiquark.) Since all of these four-quark operators couple to the scalar tetraquark state \t), we 
define the ("optimal") tetraquark interpolator J q i q i as the particular linear combination which has 
maximal overlap (0 \J q 2 q 2 \ t) with the tetraquark ground state. 




(2.1) 



m^ 2 = A(A-4) 



(2.2) 
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This definition encodes part of the tetraquark's structure into the field composition and anoma- 
lous dimension of /=2 q i since both are affected by the dominant quark couplings in The anoma- 
lous dimension y(ju), in particular, carries information on the multiquark correlations responsible 
e.g. for tetraquark binding [5]. Owing to the boundary condition (2.1) this information gets trans- 
lated into a bulk mass correction Am$ (z) for the tetraquark's dual mode 1 . (The z dependence of 
A/M5 is inherited from the dependence of 7 on the renormalization scale jU ~ 1/z.) 

Despite their rather straightforward origin, the implementation of these bulk-mass corrections 
into bottom-up duals is not without challenges. To begin with, at present QCD information on 
anomalous dimensions of hadronic interpolators is still scarce (especially in the infrared) and their 
renormalization scheme dependence raises compatibility issues. A more fundamental concern [5] 
is the naive AdS/QCD extrapolation of the multiquark physics with its pronounced N c dependence 
to N c = 3. Major benefits of the outlined mechanism, on the other hand, include that it will work 
in any of the current AdS/QCD duals, both non-dynamical (i.e. not backreacted) and dynamical 
(see e.g. Ref. [9]), and that the results are essentially model independent. In the following we will 
implement this mechanism into the "dilaton soft- wall" dual of Ref. [10] and apply it to the analysis 
of tetraquark correlations inside the lightest scalar mesons. 

3. Light scalar tetraquarks in the dilaton soft-wall dual 

The lightest scalar meson nonet (for reviews see e.g. [11]) has long been suggested to contain 
a dominant tetraquark component [1]. This at present arguably favored interpretation requires an 
exceptionally strong four-quark binding which should be large enough to reduce the tetraquark 
mass considerably below that of the lightest scalar qq mesons [12]. Hence light scalar tetraquarks 
furnish a particularly challenging testing ground for the holographic mechanism discussed above. 

To implement the latter into an explicit dynamical framework, we adopt the soft-wall dual of 
Ref. [10] which contains the AdSs metric and a quadratic dilaton field 2 <J> (z) = \ 2 z 2 ■ After Fourier 
transforming the x dependence of the tetraquark's dual modes <p and defining a "reduced" scalar 
field (q,z) = (R/z) 3 / 2 e~^ 2z2 / 2 (p (q,z), the radial bulk equation in this background turns into the 
Sturm-Liouville problem 

[-d 2 + V(z)]$(q,Z.)=q 2 $(q,z) (3.1) 

with the potential 

V(z)=(j + m|/? 2 ) i + A 2 (A V + 2) . (3.2) 

For constant the normalizable solutions of this eigenvalue problem can be found analytically 
[5]. The resulting, discrete square-masses of the 71-th radial excitations (related to A by Eq. (2.2)), 

M 2 = ql = 4X 2 (n + ^\ (3.3) 

lie on linear trajectories and indeed grow with A. 

In addition, the spectra (3.3) reveal that the ground state of the four-quark trajectory with 
Ag2 q i = 6 has a twice larger square mass M| =6 = 2M 2 q than its Aq q = 3 counterpart (as long 

'For anomalous-dimension induced bulk-mass corrections in a different context see Ref. [8], 
2 in contrast to the "metric soft wall" of Ref. [13] 
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as anomalous dimensions are ignored). Simply replacing the conventional quark-antiquark inter- 
polator by a q 2 q 2 interpolator would therefore result in four-quark states which are heavier, not 
lighter than the qq ground state. Of course, in view of our above discussion it is not surprising that 
the attraction required to bind tetraquarks is lacking. Indeed, if multiquark correlations generate 
the exceptional lightness of the tetraquark ground state, the latter should emerge from the so far 
neglected anomalous dimension of the tetraquark interpolator. 

Inclusion of this anomalous dimension, which we denote by y(z), generalizes the bulk mass 
term (2.2) to 

m 2 5 {z)R 2 = [6 + y{z)][2 + y{z)} (3.4) 

for the modes dual to the tetraquark interpolator J& q 2 with A q i q 2 = 6 + y(z). This adds the universal 
correction 

AV(z) = r(z)[7(z) + S}\ (3.5) 
z 

to the potential (3.2) with ni^R 2 = 12. Note that this correction will be completely determined as 
soon as QCD information on the RG flow of y will become available. In the meantime, we resort to 
estimating its quantitative impact on the tetraquark spectrum and to deriving an exact lower bound 
on the ground-state mass. For the former purpose, we adopt the ansatz 3 y(z) = — az 11 +bz K whose 
coefficients turn out to be tightly constrained by consistency and stability requirements [5]. 

A fundamental property of Eq. (3.5) is the lower bound AV (z) > — 16/z 2 which holds for 
any y and prevents the collapse of dual modes into the AdSs boundary. This bound is saturated by 
y = —4 and therefore yields the lower bound 

M M > M A=2 , = 2A (3.6) 

which determines the lightest tetraquark mass which Eq. (3.4) can generate. Moreover, for constant 
values of y in the range —4 < y < — 3 the tetraquark's ground-state mass remains below its qq 
counterpart. Reassuringly, the underlying binding mechanism will work in other AdS/QCD duals 
as well. Indeed, y enters the bulk dynamics exclusively through the mass term (3.4) which is model- 
independently prescribed by the AdS/CFT dictionary and yields the universal correction (3.5) in 
any AdS/QCD dual. 

Higher-lying tetraquark excitations with masses M q 2 q 2 n (for n > 0) beyond those of the corre- 
sponding quark-antiquark mesons require in addition a suitably running y(z). Our above ansatz, 
for example, can push the tetraquark masses M q 2 q 2 n beyond the M qqn from around n > 2 while si- 
multaneously generating an almost maximal ground-state binding energy. The latter drives M q 2 q 2 
from ~ 40% above (for y = 0) down to ~ 20% below the qq ground-state mass 4 M qq fl = \/6X . This 
can be seen in Fig. 1 where the corresponding spectra and potentials are plotted. It further appears 
that the higher-lying radial tetraquark excitations will be broad enough to prevent the appearance 
of supernumeral scalar states in the meson spectrum. 

3 Anomalous dimensions of a qualitatively similar z dependence (in the region of interest) are encountered in dual 
backgrounds of holographic RG-flow type [14]. 

4 Hence the phenomenological ratio Mg2„2 q/M^q ~ 0.8/1.5 is not quantitatively reproduced, perhaps due to the 
neglect of the qq interpolator's anomalous dimension. 
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Figure 1: Left panel: the tetraquark's dual-mode potential (3.2) (i) for the anomalous dimension y(z) = 
— az^ + bz K with a = 4, b = 0.05, 77 = 0.001 and K = 2 (thick line), (ii) for the four-quark interpolator with 
A = 6, y = (dotted line), (iii) for the qq interpolator with A = 3 (dashed line) and (iv) for A = 2 (thin 
line) which saturates the Breitenlohner-Freedman bound. Right panel: The ground state (« = 0) and the first 
five excitations of the tetraquark mass spectrum (i) in the potential (i) of Fig. 1 (thick bars), (ii) for the qq 
interpolator (dashed bars) and (iii) for the A = 6 interpolator (dotted bars) with 7 = 0. 

4. Summary and conclusions 

We have presented an essentially model-independent extension of the AdS/QCD framework 
which describes multiquark correlations in hadrons by means of anomalous-dimension-induced 
bulk-mass corrections. This mechanism yields a robust holographic description of a dominant 
tetraquark component in the lowest-lying scalar mesons. Moreover, it provides an explicit lower 
bound on holographic ally attainable tetraquark masses. In the dilaton soft-wall gravity dual this 
bound sits at twice the infrared scale and can be almost saturated. The resulting, exceptionally 
strong four-quark binding renders the tetraquark ground-state nonet about 20% lighter than the 
lowest-lying scalar quark-antiquark state. Higher-lying tetraquark excitations, on the other hand, 
can be pushed beyond their quark-antiquark counterparts and are thus likely to dissolve into the 
multiparticle continuum. 
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